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Abstract 


An analysis of test system misalignment is presented for the creep 
test. Sheet type rectangular 1100-0 aluminum specimens are used for discus- 
sion. It is found that the creep strain at the geometric centerline of the 
specimen is different than that at the neutral axis. However, this differ- 
ence in the creep strain decreases with time. Generally, the effect of mis- 
alignment decreases with creep time. 

Creep tests conducted with long pullrods and large initial strain level 
(high creep stress) will tend to minimize the effect of misalignment. 
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I. INTRODUCTION 

The effect of misalignment in a tensile test system has been previously 
investigated by a number of investigators. These investigators have also 
identified the sources of misalignment as poor conformance of specimen 
centerline to top and bottom grip centerlines, poor alignment of the top and 
bottom grip centerlines, and inaccurate machining of the test specimen it- 
self. A more detailed description of test system misalignment and the 
discussion concerning the need for such investigation may be found in Penny 
and Leckie [1], Christ and Swanson [2], and Wu and Rummler [3]. 

The misalignment in the creep test was previously investigated by Penny 
and Leckie [1], Hayhurst [4], and Penny et al. [5], These research workers 
generally agreed that data scatter in creep testing could be partially 
attributed to specimen alignment. Poor specimen alignment could result in 
bending stresses (or strains) superimposed on the required mean axial stress 
(or strain) . This bending stress (or strain) could be high enough to cause 
premature failure in brittle materials; in ductile materials, plastic 
deformation could occur and reduce the bending stresses (or strains). 

Analyses of misalignment were presented in Ref. [1 and 4] for the 
creep test. The results show that the creep curves "scatter" due to eccen- 
tric loading. More significantly, these analyses predicted that the 
"scatter" would increase with the creep time. 

The result of the present analysis does not support all of the above 
assertions. While the numerical results show that the errors induced by 
the test system misalignment are significant in the creep test, these errors 
will diminish with time. Thus, the most severe effect of misalignment 
occurs at the beginning of the creep test. 
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The analysis Is performed on rectangular specimens of the sheet type. 
Specimens of 1100-0 aluminum Lave been chosen for investigation. The 
reason for using this material as an example for discussion is that the con- 
stitutive equation employed in the present analysis has been shown to 
describe both the dynamic stress-strain behavior and creep for this material 
[6] and the material constants have already been determined. Theoretically 
speaking, the method presented herein would also be applicable to other 
materials. 

Only Uie symmetric case of misalignment is studied due to the finding 
of Ref. [3] that this is the most severe case of all the cases investigated. 

The endochronic constitutive equations developed by Wu and Chen [6] for 
creep are used in the analysis. However, it is believed that the results 
obtained in this investigation concerning the effect of misalignment is 
independent of the constitutive model used. Any consitutive creep equation 
which can correctly describe the creep behavior may be used for the mis- 
alignment analysis. 

The present analysis is the third part of a series of investigation 
related to test system misalignment. The first part was presented in Wu 
and Rummler [3] which reported a misalignment analysis for the tension test 
under static loading condition. The second part was given in Wu et al. [7] 
and reported an analysis of the tension test under dynamic loading condition. 
In particular, the effect of strain-rate was investigated in Ref. [7]. The 
results of this work also serve as the initial input to the present calcula- 
tion. 
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II. CONSTITUTIVE EQUATION FOR CREEP 

The present analysis utilizes a constitutive equation which was shown 
by Wu and Chen [6] to describe the transient creep and creep recovery. In 
this approach, the deformation history of the specimen prior to the onset 
of creep is important and significantly affects the transient creep behavior. 
Therefore, the whole creep test is viewed to have two stages. The first 
stage is the loading stage and the second stage is the creep stage (see 
Fig. 1). 

The loading stage is treated as a constant strain-rate tension test 
(see Ref. [6] for details). The numerical results reported in Ref. [7] are 
precisely the numerical data for this stage of the test. The creep stage 
is described by the following equations (see Ref. [6]): 
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In the above equations, o q , 6, n and 8 R are material constants de- 
termined by a set of constant strain- rate stress-strain curves; o* is the 

creep stress which is a constant; 9 is the creep strain; k and z are 
r 00 

values determined at the end of the loading stage; and 6 depends on a* as 

c 

a parameter. 

Rearranging equation (3), the creep rate is expressed by 
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( 10 ) 
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Thus, equations (1), (2) and (10) form the basic equations describing the 
transient creep. 

To compute the creep curve, a step-by-step numerical scheme should be 

employed. Generally, the current value of k^(0) is computed from equation 

(1) by use of the value of z at the previous time step. The current value 

of k (0) is then kept constant when the increments of the creep strain and 
c 

the value of z are being computed from equations (10) and (2). 
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III. ANALYSIS OF MISALIGNMENT 

As mentioned in the previous section, the creep test under investigation 
is a constant strain-rate loading test followed by a constant stress creep 
test. In the loading stage the strain-rate is constant at the neutral axis 
(N.A.) of the specimen. In the creep stage, the stress along the N.A. 
remains constant, but it changes with time at other parts of the specimen 
due to the effect of misalignment. 

In fact, phenomena such as creep recovery and stress relaxation may 
indeed occur at the part of cross-section outside of the N.A. However, the 
stress variation for a given fiber at each time step of computation is so 
small that the effect of creep recovery may be ignored compared with the 
creep strain resulted from the constant stress. This assumption will be 
justified by the numerical results. 

The assumption of plane cross-sections remaining plane during deforma- 
tion is made in the analysis. The coordinate system is same as that used 
in Ref. [3 and 7] and shown in Fig. 2a. The essential features of the sym- 
metric case of misalignment is also presented in Fig. 2b, The discretized 
numerical model is identical to that of Ref. [7] and is reproduced in Fig. 3. 
Due to symmetry, only half of the specimen is considered, which is then 
divided into I and J sections along the longitudinal and the lateral 
direction, respectively. 

At a particular time t in creep, the planes remaining plane assumption 
requires that 
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where x and y are coordinates shown in Fig. 2(a); e is the strain on the 
N.A.,; p(x) is the radius of curvature of the specimen; e and £ are the top 
and bottom fibers, respectively; h is the width of the specimen; and d(x) 
specifies the location of the N.A. and is measured from the bottom fiber. 

The superscript t indicates that the quantities are referred to the time 
step t. 

The procedure for solution is such that the quantities at the N.A., 
represented by the barred quantities, should be assigned first. This may 
be done by use of equations (1), (2), and (10). Therefore, 

e t - e R exp {(i - 6 c > (14) 
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The creep behavior at the N.A. is same as that of a perfectly aligned 
specimen. Quantities at other fibers are then determined by means of the 
misalignment analysis given in the remaining part of this section. 
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At a given time t, balance of load is 

P - f a t dA - ^ l aj ; i - 1 1+1 (17) 

A J j-1 iJ 

in which P is the force, b is the thickness, A is the cross-sectional area, 
and h/J is the width of the element under consideration at the i*" section. 
The balance of moment is given by 
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where equation (11) was used. Furthermore, Fig. 2(b) shows that the moment 
at the i 1 "^ plane may be related to the eccentricity 6 by the equation 
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deflection, L is the length of the pul 1 rod, and 
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Equations (18) and (19) combine to yield 
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In addition, the curvature relation is 
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( 22 ) 


It is then easily shown that the following equations may be obtained 
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where E is Young's modulus. 

On the other hand, based on the discrete model, the creep constitutive 
equation should be used to correlate the stress and strain for the element 
(i,j). For this element, equations (1), (2) and (10) are written as 
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Note that the "hatted" quantities are those related to constitutive creep 
equations. These quantities should be matched by Iteration with the "un- 
hatted" quantities which appeared in the field equations. 

An assumption is made in arriving at these equations. It is assumed 
that the creep strain at a generic point in the specimen may be calculated 
by maintaining the same stress at this point throughout the creep stage. 
Thip, in effect, neglects the effort of creep recovery and stress relaxation 
at fibers other than the N.A.. 

The matching of the above mentioned quantit-.a.*. are represented by the 
equation 


ej - ej ; i - 1, .... 1+1 and j - 1,..., J (32) 

This equation also represents the Input by the constitutive creep equation 
into the formulation of the misalignment problem. Hence, the set of non- 
linear equations (17), <„20) , (21), (22), (27), and (32) represents the basic 
equations of the misalignment analysis. The unknown variables are c* e*g. 
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£ t t t 

d g , and and, in the computation, these equations should be sup- 
plemented by equations (12), (13), (23), (24), (25) and (26). In solving 
the set of nonlinear equations, the right hand side of equation (32) is 
considered known. 

The set of nonlinear equations may be solved with the following 
boundary conditions: 


— - « 0 at x - 0 
dx 


(33) 


and 


u fc = 0 at x * Jt/2 (34) 

provided that the initial values (at t = 0) of all the variables which cor- 
respond to a particular value of 9 q are supplied from the results of Ref. 
[7]. 

An algorithm for the numerical solution of the above set of equations 
is presented in the Appendix. The theoretical results are discussed in the 
following section. 
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IV. RESULTS AND DISCUSSION 

The sheet type rectangular specimens are used in this calculation. The 
geometrical data for the specimens are 

specimen length (£) « 63.5 mm 
specimen width (h) = 12.7 mm 

specimen thickness (b) « 3.2 mm 
eccentricity (6) * 0.05h 

The influence of the pull rod length (L) is investigated by assigning two 

different ratios of L/£, i.e., L/£ = 2.0 and 1.25. 

The influence of the preloading history on misalignment is investigated 

through the following consideration. The specimen is first subjected to a 

• * -5 -1 

loading process of constant strain rate 0 “0^“1.3OxlO S at the N.A. 

The loading continues until a prescribed plastic strain 0 Q is reached at the 
N.A. The load is then kept constant and creep is subsequently computed. 

By considering different values of 0 , the influence of prehistory may be 
studied. 

In this calculation, three cases are considered, which are L/£ = 2, 

0 * 0.3%; L/£ - 2, 3 = 0.5%; and L/fc = 1.25, 0 » 0.5%. Results related 

o o o 

to the influence of the pull rod length and the prehistory will be discussed 
later in the text. The initial data at the preloading stage, such as the 
initial values of misalignment errors and the location of the N.A. prior to 
the creep stage are supplied by the misalignment program for the dynamic 
tension test reported in Ref. [7], 
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In order to reduce computer time, only four sections have been cut along 
both the longitudinal and the lateral directions of the specimen, l.e., 

I * A and J ■ A. Thus, the unknown variables have been reduced to 55 in 
number and the set of nonlinear equations discussed in the previous section 
consists now of 55 equations. Although the number of elements in the dis- 
crete numerical model is at its minimum, the model does provide consistent 
numerical output under a variety of creep conditions. Since the area of 
interest for the misalignment analysis lies in the midspan of the specimen 
(the measuring devices are usually attached to the midspan) , most of the 
results to be discussed are related to this area. 

The material constants used in the calculation have been determined in 
Ref. [6]. They are 

o° » 13.79 MPa (2000 psi) 

o° « 27.92 MPa (A050 psi) 

E l - 6.596 x 10 3 MPa (9.566 x 10 5 psi) 

n = 105 
B » A. AAA 

The parameter of equation (3) has been found in Ref. [6] to depend on 
the creep stress. This dependence may be represented by the following 
linear equation (Fig. A) 

8 . - 3.6A3 x 10' 2 - 0.261 x lO -1 (a* - 0.508) 

cij ij 


for o* > 0.A93 MPa 


(35) 
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where the creep stress is measured In MPa. For stress smaller than 0.493 
MPa (or 3400 psi), the relationship between $ c and o* Is non-linear. How- 
ever, in all the cases considered in this calculation, the creep stress is 
in the range of equation (35) and the linear equation is used in the computa- 
tion. 

Figure 5 shows the normalized stress a*^/a* plotted versus time for the 
top and bottom fibers in the three cases considered. It is seen that in 
all cases the variation in stress is minimal. The change in stress is less 
than 1% in the first hour of the creep test and is less than 2% in 8 hours. 
Therefore, the assumption made earlier that the effect of creep recovery 
and stress relaxation are negligible is a reasonable one. Figure 5 further 
shows that it is unlikely that a uniform stress distribution in the specimen 
can be reached before rupture. 

It may also be observed from Fig. 5 that, for a fixed L/SL ratio, the 

deviation of the normalized stress from unity at any time t is greater 

for smaller 0^. The implication is that greater bending moment is obtained 

for smaller 0 . This observation is related to the inherent error range 
o 

for the extreme fibers at the pre-loading stage reported in Ref. [7J. A 
conclusion of Ref. [7] is that the error range for the extreme fibers is 
always decreasing with the increasing 0 . Figure 5 also indicates that this 
error range in stress for a fixed 0 q decreases with the creep time. 

A second observation from Figure 5 is related to the length of the pull 
rod. The longer the pull rod is the smaller is the error range in stress. 
This observation is closely related to the magnitude of the end moment which 
will be further discussed. 


% 
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In Figure 6, the creep strains versus time for the top and bottom 
fibers at the midspan of the specimen are shown and compared with the creep 
curves at the N.A. Although the N.A. shifts with time, it is convenient 
lo a.';-.; the creep response at the N.A. as a reference and compare those at 
otht • fibers with it. It is seen from Fig. 6 that for a given 6 o> a shorter 
pull rod will result in a greater variation in the magnitude of strain 
between the two extreme fibers. Furthermore, for a specific pull rod 
lengt l, the creep rate is smaller for smaller 8 q . 

lit appears in Fig. 6 that for all the three cases considered the creep 

rate of each fiber were almost identical to that of the N.A, and all creep 

curves, for each case were parallel to each other. A closer look at these 

curves will show that this is not true. Indeed this is made clearer by 

Fi 0 s. 7 and 8, in which the differences in strains between the extreme 

fibers and between the bottom fiber and the N.A. are plotted against time 

t. These curves show two different trends of variation. For 9 = 0.5%, 

o 

the differences : l strains reach their peak values within one hour and 
decrease with time thereafter. On the other hand, in the case of 0 q = 0.3%, 
the strain differences increase with the creep time up to 8 hours. Due to 
the limitation in computer funds, the computation is terminated at 8 hours 
creep time. But i" is conjectured that these curves, would reach their 
peaks at a late time and show a downward trend thereafter. Thus, the 
peaks wou. 1 . occur early for a large 0 q . 

F gure 9 shows the difference in strains between the geometrical center- 
litii and the N.A. at the midspan. It is seen that these differences in 
strain are always decreasing with time for all the cases considered. However, 
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the computation time Is not long enough to see If the differences will 
eventually vanish with time. Nevertheless, It may be concluded, that for 
the same pull rod length, the smaller the prestrain Q q Is the longer Is the 
creep time required to close the difference between e c (the strain at the 
centerline) and e (the strain at the N.A.). 

Perhaps the most significant of all the results of this paper is that 
presented in Fig. 9. The implication of this result is that the difference 
between e c and e will decrease with time in spite of the test system mis- 
alignment. That is, the creep strain measured by an extensometer, which is 
also the strain at the geometric centerline e^, will become closer and closer 
to the creep strain of a perfectly aligned test system at large time. 
Therefore, it may be concluded that the effect of misalignment in a creep 
test is significant only in the primary creep range. In a long time creep 
test, the effect of misalignment is less important, because the plastic de- 
formation will even out the non-uniform distribution of strain with time. 

This finding is not in agreement with those of Penny and Leckie [1] and 
Hayhurst [4J. The analyses of Ref. [1] and [4] predicted that the eccen- 
tricty cf the test system would cause the creep curve to deviate increasingly 
with time from the creep curve of the perfectly aligned test system. 

The error of misalignment at the geometric centerline (or center error) 
is defined by = (e c - e)/e. The center errors at the midspan for the 
three cases under consideration are shown in Fig. 10. It is seen that the 
center error decreases with time for all cases. The errors are within 2% 
at the onset of creep and become less than 0.4% after 8 hours of creep. 

The same trend is found in the curves of Fig. 11 for the errors of top and 
bottom fibers. These errors are within 30% at the onset of creep and are 
about 8% to 12% after 8 hours of creep. 
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The shifting of the N.A. and the end moment distribution are closely 
related by equation (19) and are plotted in Figs. 12 and 13 against time. 

All curves are of the decreasing type with time. However, a conclusion 
has not been reached concerning the questions of whether the N.A. would 
eventually coincide with the geometric centerline and whether the end 
moment would eventually be reduced to zero at large time. More computations 
are needed in order to answer the above questions. 

Finally, the deflection curves of various creep times are plotted in 
Fig. 14. For the cases of L/i*2, 0 q “ 0.5% and L/i, * 1.25, 0 q “ 0.5%, 
the deflection curves reach their maximum position within 2 hours and then 
decrease their magnitudes thereafter. For the case of L/i ■ 2 and 8 q * 0,3%, 
the magnitude of the deflection curve increases with time for the 8-hour 
creep time period. It is conjectured that the deflection would reach its 
maximum position at a later time and the magnitude would then decrease with 


time. 
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V. CONCLUSIONS 

The formulation of a misaligned specimen subjected to the creep test 

has been presented. The investigation is restricted to the symmetric case of 

test system misalignment. The specimens used in the calculation are the 

sheet type rectangular specimens of 1100-0 aluminum. Three cases have been 

investigated in detail. They are the case of L/i « 2 and 0 q - 0.3%, the 

case of L/fc * 2 and 0 » 0.5%, and the case of L /l * 1.25 and 0 ■ 0.5%. 

o o 

The following results have been obtained from this investigation: 

(1) During the creep test, the stress is quite constant for all fibers. 
The total change in stress for the 8-hour creep period for the top and bottom 
fibers is within 2% for all cases considered. 

(2) The difference between the creep strains at the geometric center- 
line and the N.A. decreases with the creep time in spite of the test system 
misalignment. 

(3) The center error and the errors at the extreme fibers are always 
decreasing during creep. 

(4) For the same pull rod length but with different amounts of pre- 
strain, the smaller the prestrain is (in the plastic range) the greater is the 
misalignment error at the onset of creep and the longer creep time it takes 
fjr the misalignment effect to vanish. For the same prestrain but with 
different pull rod lengths, the error is greater for the shorter pull rod 
than for the longer. 

(5) In order to effectively reduce the misalignment error in a creep 
test, long pull rod should be uc'id and the test should be conducted with 

a large prestrain. 
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APPENDIX 

Algorithm for Numerical Solution of the 
Misalignment Problem in Creep 


In this calculation, the specimen is discretized into 4 elements along 
both the longitudinal (1*4) and the lateral (J*4) directions. Algebraically, 
there are 55 nonlinear equations with 55 unknown variables. They are: 


Equations: eq. (17) - 5, eq. (20) - 1, eq. (21) - 5, eq. (22) - 4, 

eq. (27) - 20, eq. (32) - 20. 

Variables: o^ - 20, e* T - 5, e^ B - 5, z^ - 20, u* - 4, d* - 1. 

The 55 variables are denoted by X(55). The computation is designed to 

achieve the following three tasks: (a) To compute the strain £ ^ in the 

element (i,j) from the constitutive creep equations, (b) To match the strain 

with the strain and solve the 55 nonlinear equations, (c) To 

iterate the step (b) until and * o*^ . 

Step (a) includes setting up m (approximately .00) creep stress levels 

between the extremum stresses in the specimen. These are generally the 

stresses at the top and bottom fibers of the grip end. The m stress levels 

should include the creep stress o* at the N.A. The initial values (at t * 0) 

7 ° k°, e° for each stress level are found by iteration. The m creep strains 

m m m 

for each time step are then found from equations (28), (29), (30), (31) and 
(35) with the subscripts ij replaced by m in these equations. 

In Step (b) , the initial guess of the variables is modified by 


X 


t 


x £ -“ (s*/i t_4t ) 


(36) 


I 1 


I 

! 

I* 


i 


* 
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The 20 stresses X(l) ^ X(20) are then matched with the m stress levels found 
in Step (a). Linear interpolation may be used between stress levels. Thus, 

A £ 

the creep strains may be obtained. 

* t 

Putting the values of found in Step (b) into the right hand side of 
equation (32), the set of 55 equations is then solved numerically. However, 
the stresses found may not agree with the 20 stress levels which give 
rise to Thus, in Step (c) an iteration procedure is employed to 

achieve the goal that ■ o*j when . 

The above procedure is repeated for t + At. 
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FIGURE CAPTIONS 
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7 Strain Difference Between Extreme Fibers 

8 Strain Difference Between Bottom Fiber and N.A. 

9 Strain Difference Between Geometric Centerline and N,A. 

10 Center Error vs. Creep Time 

11 Errors at Top and Bottom Fibers 

12 Shifting of N.A. 

13 End Moment 

14 Deflection Curves 
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